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A numerical calculation of properties of finite absorbers and of intergalactic 
medium based on photoionization equilibrium is performed to confront alter- 
native UV sources in addition to quasars. 

It is seen that a spectrum including a large peak around the HI ionization 
energy due to decaying neutrinos is too soft in the region up to the Hel edge 
i to explain the relatively small observed ratio of neutral He and H densities in 

_ Lyman-limit systems if their size is of the kpc order. The recently proposed 

■ decrease of the contribution from unstable neutrinos solves this problem but 

. tends to spoil the consistence between IGM and Lyman-a clouds, which requires 



> 



ON 



a large ratio of fluxes for the HI and Hell ionization frequencies, unless there 



• is a very fast decline of quasars above z = 3. 

On the other hand, the addition of stars to quasars may produce a spectrum 
O I sufficiently hard between HI and Hel and thereafter soft up to Hell to allow 

' a reasonable agreement of the properties of denser absorbers with those of 

■ IGM. This model seems to favour cold dark matter with additional cosmological 

constant. 

^ ' Key words: dark matter - galaxies: intergalactic medium - quasars: absorption lines 

1. Introduction 

The reionization of the Universe is a subject of deep experimental and theoret- 
ical study for its implication on galaxy formation and on the consistence of small 
scale anisotropies with cosmological models. The different possible sources of UV 
background may give some information on the identity of dark matter. 

It has been seen from the bounds on Gunn-Peterson (GP) effect (Gunn & Peterson 
1965) for neutral H and single-ionized He that the frequency UV spectrum seems to be 
softer than that due to quasars (QSO) (Madau 1992), requiring therefore additional 
sources which might be stars. 

Another suggested possibility has been that decaying neutrinos (Sciama 1990a), 
in addition to QSOs, may ionize HI in the intergalactic medium (IGM), clouds and 
other systems. This decaying dark matter (DDM) would be hot (HDM) and might 
be in difficulty to explain the formation of structures unless cosmic strings are the 
seed of primeval fluctuations (Zanchin et al. 1996). Since this is a delicate subject, 
it is interesting to see whether the sole properties of reionization are able to support 
this hypothesis or not. 
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The observations whose vahdity we will assume for our analysis are the ratio 
of densities of neutral He and H in Lyman-hmit systems (LLS) (Reimers & Vogel 
1993), the bound of GP effect for HI (Giallongo et al. 1994) considered as due to 
photoionization and the estimation of the same effect for Hell (Jakobsen et al. 1994; 
Tytler et al. 1995; Davidsen et al. 1996). 

We will calculate the He ionization fractions for absorbers of given density due 
to UV fluxes corresponding to the different alternatives to inspect their agreement 
with the observations of LLS. Additionally, from the expressions of GP effect for HI 
and Hell we will obtain a bound for the ionization fraction of the latter in IGM to 
compare its consistence with that of Lyman-a clouds (LC) for the various possibilities 
of flux. 

In Section 2 we will describe the possible sources of UV radiation whose difference 
is roughly that QSOs ionize HI, Hel and Hell, stars just HI and Hel and decaying 
neutrinos only HI. The contribution of stars will be fixed by recent determinations 
of the proximity effect (Giallongo et al. 1996) and that of DDM requiring that it is 
capable to ionize alone the HI and NI of the Milky Way (Sciama 1990b) or that the 
decay photons have not enough energy to ionize NI (Sciama 1995). 

In Section 3 the numerical calculation will be presented for absorbers of different 
density with the alternative UV fiuxes to determine the ionization of He and other 
properties and see their compatibility with observations in LLS. 

Section 4 will be devoted to compare GP effect for HI and Hell with the different 
models estimating the fraction of the latter in homogeneous IGM to evaluate through 
approximate formulae if it is consistent with the numerically calculated properties of 
LC. 

The conclusions will be given in Section 5 indicating the possible relation of the 
UV ionizing flux with different cosmological models. 

2. Sources of UV background 

Recent determinations of bounds for GP effect for HI and estimations for Hell 
give, under the assumption of photoionization to explain the ratio of the optical 
depths 

r§^n/r§f = 0.45 Jm/JHeii , (1) 

a ratio of the effective UV fluxes at the corresponding ionization frequencies of at 
least 

Sl = Jhi/JhcII = 100 (2) 

for z — 3.3. Even if the GP effect is masked by line blanketing, the above ratio is 
maintained rather high with the estimation Sl > 40 (Madau & Meiksin 1994) or 
Sl > 65 (Sethi 1995). On the other hand, from metallicity abundance it is found 
that 5l - 70 at z ~ 3.2 (Songaila et al. 1995) and even Sl > 100 for z ^ 3.5 - 3.8 
(Savagho et al. 1996). 

This ratio Sl seems to exceed that due to QSOs which was evaluated as ~ 30 
(Madau & Meiksin 1994; Haardt & Madau 1996) in agreement with the declining 
quasar population for 2; > 3 (Pei 1995). It appears therefore necessary to fill the 
difference with another UV source. 



9 



If one adds stars of primeval galaxies (Miralda-Escude & Ostriker 1990) they may 
ionize HI and Hel with a spectrum similar to that of QSOs but which afterwards drops 
abruptly so that Hell is ionized only by the latter source. If one fits the resulting flux 
of QSOs in the range between HI and Hel as J ~ v~°' with o; < 1, it follows that for 
the total flux 



Another possibility is the addition to QSOs of DDM thought (Sciama 1993) as 
neutrinos of mass around 30 eV which could correspond to close the universe with 
HDM. This model is alternative to the first since HDM delays galaxy formation so 
that stars would be excluded. According to the details of the decaying neutrinos the 
relation of Eq. 3 may be altered. 

Our analysis will be based on considering three models for UV fluxes: that due 
to QSOs alone and the alternatives of adding cither stars or decaying neutrinos. We 
will consider them at z — 2 and 4 because the former redshift corresponds to the 
observations of LLS, and the latter to the range of GP estimations for HI and Hell 
(see Fig. 1). 

Using for the effective UV flux the normalization 
J = J_22 X 10~^^ erg s~^ cm"^ Hz~^ sr~^, the evaluation of the QSO contribution 
has increased doubling the earher values (Madau 1992) to reach at IRy a maximum 
of J_22 ~ 3 for ^ > 2 when dust obscuration is taken into account (Haardt & Madau 
1995) and diminishing for larger redshifts so that we estimate J_22 ~ 2 at 2; ~ 4. 
According to the discussed values of Sl due to QSOs, and considering that this ratio 
increases with z, one evaluates as maximum contribution at iRy J_22 ~ 0.05 at 
z — A and shghtly above 0.1 at 2; = 2, which will coincide with the total flux at this 
frequency in all the alternative models. We will denote this scheme as the strong 
QSO model (sQSO). 

Stars give a slightly softer spectrum in the range HI - Hel and we normalize 
the total fiux in the QSO+star model at IRy as J_22 ~ 5 to agree with recent 
determinations from proximity effect J_22 = 5 ± 1 at IRy which show no evolution 
in the range 2 < z < 4 (Giallongo et al. 1996). In this way it will be noted that the 
ratio of Eq. 2 is satisfied for 2; ~ 4 and that of Eq. 3 reasonably fulfilled in the range 
2 < ^ < 4 . To explore the situation Sl = 200 at ^ ~ 4 we will also consider the weak 
QSO contribution (wQSO) corresponding to take 50% of the above quoted values, 
always keeping the normalization J_22 = 5 at IRy for the total QSO+star flux. 

It must be remarked that the fluxes denoted as Jhi, JhsI ^-nd JhcII are the average 
over frequency weighted with the corresponding ionization cross section. In the QSO 
and QSO+star cases they differ only slightly from the flux J,, at u = IRy etc, because 
of the fast decrease of the cross section. E.g. 



On the other hand for the QSO+DDM model the only difference from QSOs alone 
is a large peak around IRy so that this maximum of will be much greater than 
the averaged Jhi- From the requirement (Sciama 1990b) that decaying neutrinos 
are able to ionize the H of the Milky Way their lifetime must be r,^ ~ 2 x 10^^ s. 



Jhi/ JhcI < 2 . 



(3) 




Imyi'^Hi/i^) du 



(4) 



Figure 1: The UV ionizing source as function of the frequency (in Rydberg) for two 
different redshifts. The continuum hne is for QSO sources only, the dashed hne is 
the same with the addition of stars, the dotted hne is with the addition of decaying 
dark matter (DDM) photons. The upper two plots represent the UV flux in case 
of strong contribution by the QSOs (sQSO), and in case of the QSO+DDM model, 
of strong contribution by the DDM (sQSO+sDDM). The lower two panels are for 
weak QSO contribution (wQSO) and, in the case of the DDM model, for strong QSO 
contribution, but weak DDM contribution (sQSO+wDDM). 
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Correspondingly, the intergalactic flux expressed in number of photons per cm^ and 
s at z ~ would be 



Ho 13.6 ^ ' 

with photon energy = (13.6 + e) eV, and where Ui, is density of neutrinos and Hq 
present Hubble constant. Assuming (Sciama 1993) that decay photons are able to 
ionize NI, the neutrino mass must be rrii, ~ 29 eV and consequently F ~ 6 x 10^ which 
corresponds to J_22 ~ 24. Since according to this model the flux due to neutrinos 
increases with rcdshift as (1 + z)^^'^, we would have J_22 ~ 120 at z = 2 . 

The ionization cross-section has the frequency behaviour a ~ so that for a 
peak of width e additional to the spectrum due to QSOs one has, according to Eq. 4, 

T _ tDdm 3 e Qso 3 
^^^-■^ 13:6 + "^^^^ 3T^ • 
For the above values, using strong QSO contribution, one would obtain 
Jhi -22 ~ 27 at z = 2 which is extremely large. Considering a reduction of the flux 
due to absorption by clouds, one may estimate at IRy J -22 ~ 80 for z — 2, case which 
we will denote as strong DDM model (sDDM). As a consequence, the averaged flux 
for QSO+DDM will be Jhi -22 ~ 18 and we may establish a typical ratio (Sciama 
1994) at z = 2 for this standard QSO+DDM model 

Jhi/JhcI ~ 8 . (7) 

For z — A the same strong DDM contribution will be J_22 ~ 150 which gives an 
averaged flux Jhi -22 ~ 30 for sQSO+DDM and therefore a ratio Sl much larger than 
that of Eq. 2. It is seen that the flux due to neutrinos is higher than the large initial 
observations of proximity effect at 2; ~ 2.5 (Bajtlik et al. 1988) and even beyond 
the upper quoted values (Bechtold 1994). Other determinations give smaller fluxes 
which either would be consistent with QSOs alone (WiUiger et al. 1994; Lu et al. 
1996) or correspond to the value we have taken to normalize the QSO+star model 
(Giallongo et al. 1996). Therefore for the described QSO+DDM case a non-standard 
interpretation of the proximity effect must be invoked (Sciama 1991). 

But recently a modified decaying neutrino has been proposed (Sciama 1995) due 
to an estimation of the metagalactic flux F at 2; ~ four times smaller than that 
given by ~ 29 eV. Therefore, keeping the value of the lifetime, from Eq. 5 e must 
not exceed 0.2 and the decay photons cannot ionize NI. The maximum fluxes due to 
this weak DDM (wDDM) would be J_22 ~ 30 at 2; = 2 and J_22 ~ 67 at 2; = 4 which 
would be roughly compatible with the highest estimations from proximity effect. Due 
to the fact that the peak is now smaller and narrower, from Eq. 6 the averaged flux 
for the modified DDM+sQSO model turns out to be Jhi -22 ~ 4.4 at 2; = 2 and 
Jhi -22 ~ 5 at 2; = 4 . Therefore the bound (3) is almost respected and that of 
Eq. 2 exactly satisfied. We note that we have performed a more accurate relation 
between neutrino mass and flux of decay photons than that done in Miralda-Escude 
& Ostriker 1992. 

Since the standard QSO+sDDM model satisfies the Eq. 7 whereas the QSO+star 
alternative as well as QSOs alone and the modified QSO+wDDM model correspond to 
a ratio (3), one may test their consistence using them in the ratio of neutral densities 



which follows from photoionization equilibrium for HI and Hel (Miralda-Escude & 
Ostriker 1992; Sciama 1994) 



= 0.044- — XHeii ■ (8) 

nni JhcI 

XhcII is the fraction of single-ionized He and the numerical coefficient is independent 
of temperature and density of the system. Moreover, information from GP effect 
which involves HI and Hell does not support QSOs as the only UV source and may 
distinguish between QSO+star model and the modified QSO+wDDM possibility. 

3. Lyman-limit systems and Lyman-o; clouds 

For LLS a ratio UHei/nHi ~ 1/30 has been observed (Reimers & Vogel 1993) for 
z ~ 2. If one takes the standard QSO+sDDM alternative of Eq. 7, Eq. 5 would 
require that for 2; ^ 2 XhcII — 0.1 the rest being almost all Helll. On the contrary 
if one chooses either the fluxes of QSOs alone or those of QSO+star or modified 
QSO+wDDM models, it follows from Eq. 3 that XhcII > 0.5. 

To obtain the fractions of ionized He, we have performed a numerical calculation 
for absorbers using the standard photoionization code CLOUDY (Ferland 1991). 

As input for the LLS we have taken the HI column density 
Nhi ~ 10^^ cm~^, the metallicity as 1/100 of the solar value and the H density 
uh ~ 10~^ cm~^. This last figure corresponds to the lower bound (Steidel 1990; 
Lanzetta 1991) consistent with sizes smaller than 15 kpc. 

For the LC we have fixed Nhj ~ 10^^ cm~^, uh ~ 10~^ cm~^ and metaUic- 
ity 1/1000 of the solar value, which correspond to their commonly accepted sizes 
(Meiksin & Madau 1993). Lower limits to the transverse size of the LC observed in 
the gravitational lensed spectra of HE1104-1805 (Smette et al. 1995) is of the order 
of 50 kpc at 2a level for spheroidal clouds at z ~ 2, where h — Hq/(100 km 
s"^ Mpc~^) . Similar results have been found (Dinshaw et al. 1995) from the cross 
correlation of absorption lines in the close quasar pair 1343+2640A and B deriving 
sizes larger than 40 kpc at z = 1.8. 

We have assumed that the UV background radiation is the only ionizing source 
and we used for the flux the different shapes and intensities according to Fig. 1. 
Results are summarized in Tables 1 and 2. 

Considering LLS, it turns out that at 2; = 2 XhcII > 0.5 for the three models. In 
particular for the modifled decaying neutrino in addition to the sQSO contribution, 
the values of XHeii turn out to be similar to those of the sQSO+star model. Therefore 
XHeii ^ 0.1 is excluded with all the models for the chosen value of the density uh- 

One must remark that the observation of nHei/''^Hi in LLS was accompanied by 
that of ionized states of C, N and O (Reimers & Vogel 1993) which would be consistent 
with a large He ionization, i.e. XhcII ~ 0.1. This led to the explanation of the whole 
set of observations with the standard decaying neutrino (Sciama 1994) using the large 
ratio (7) of Jhi/ JhcI- But to obtain xhcII ~ 0.1 it is necessary to decrease the density 
Uh in such a way that the sizes of LLS would be larger than 100 kpc which is an 
unattractive possibility. This has suggested that the heavy elements are collisionally 
ionized in hot regions (Giroux et al. 1994) different from those where H and He are 
photoionized. This is the scenario adopted for our analysis. 

Therefore we may state that QSO+star model, QSOs alone and also the modifled 



Table 1: Fraction of Hell and Helll for different UV background models (as shown 
in Fig. 1) in the LLS. 



Model 




z = 2 






Z = 4: 






J{lRy) 


-22 XHell 


XHelll 




-22 XHell 


XHelll 


sQSO 


3 


0.60 


0.40 


2 


0.78 


0.22 


wQSO 


1.6 


0.70 


0.29 


1 


0.86 


0.14 


sQSO+star 


5 


0.67 


0.33 


5 


0.84 


0.15 


wQSO+star 


5 


0.79 


0.20 


5 


0.91 


0.09 


sQSO+sDDM 


80 


0.79 


0.21 


150 


0.94 


0.06 


sQSO+wDDM 


33 


0.67 


0.33 


69 


0.87 


0.12 



Tabic 2: Fraction of Hell and HeHI for different UV background models (as shown 
in Fig. 1) in the LC. 



Model 




z = 2 






z = A 






J(l%). 


-22 XHell 


XHelll 


J{lRy). 


-22 XHell 


XHelll 


sQSO 


3 


0.004 


0.996 


2 


0.013 


0.987 


wQSO 


1.6 


0.009 


0.991 


1 


0.029 


0.971 


sQSO-|-star 


5 


0.004 


0.996 


5 


0.013 


0.987 


wQSO+star 


5 


0.008 


0.992 


5 


0.026 


0.974 


sQSO+sDDM 


80 


0.004 


0.996 


150 


0.015 


0.985 


sQSO+wDDM 


33 


0.004 


0.996 


69 


0.014 


0.986 
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sQSO+wDDM alternative are consistent with the observation of nHei/nni aX z = 2. 
Prom the tables we may see that the CLOUDY calculation predicts a smaller He 
ionization at 2; = 4 as it is generally expected. 

For further considerations it is convenient to use approximate expressions assum- 
ing photoionization equilibrium for HI and Hell in a thin absorber of size R expressed 
in kpc and column density given in cm~^, which allow to write (Miralda-Escude & 
Ostriker 1992) 

( Nhi J hi, -22 \ 2.7 X 10-1° 

XHell - r-^ —J , (9) 

where the numerical figure includes the square root of the recombination coefficient 

whose temperature dependence is (Spitzer 1978) a{T) ~ 7"-3/4 g^j^j^ defining 

T — T4 10'' K, T4 2 has been taken. It is possible to check the densities for 

which the approximate expressions are valid so that also IGM can be included in the 

comparison. 

For a thicker absorber like LLS the left-hand side of Eq. 9 must be replaced in 
principle by XHeii/i^ — XHeii) since the right-hand side gives nueii /nHeiii- But 
in this way the obtained XHeii will be a lower bound because self-shielding for He 
should be already included for LLS, diminishing the value of JReii- Another source 
of approximation of Eq. 9 is the fact that instead of the detailed fiux spectrum, the 
averaged values are included what is particularly relevant in Jhj for the QSO-I-DDM 
case, as indicated in the previous Section. 

Comparing with our more accurate results of Table 1 for LLS, it is easily seen 
that the approximate formula gives typically a value of XHeii lower in 0.1 for QSOs 
and a slightly larger difference for QSO-|-star and QSO-I-DDM models. In fact for 
z = 2 CLOUDY gives e.g. for sQSO, sQSO+star, sQSO+sDDM and sQSO+wDDM 
models (Table 1) sizes R ~ 1.4, 2.4, 5.2 and 2.2 kpc respectively, with a temperature 
T4 ~ 2. The difference can be understood if in the approximate formula JhcII is 
reduced around 40% in the QSO case because of self-shielding, around 50% in the 
QSO-Fstar and QSO+wDDM models and even more than 60%. for QSO-FsDDM due 
to the increasing size of the absorber which results from CLOUDY. 

For thin absorbers like LC, one expects that both CLOUDY and the approxi- 
mate Eq. 9 give similar results because the absorber becomes thick to Jneii only at 
Nhi ~ 10^^ cm-2 (Miralda-Escude & Ostriker 1990). In fact CLOUDY gives for 
LC a ionization larger than for LLS with a typical XhcII > 0.01 at z — A (Table 
2). One may check that Eq. 9 gives the same values of Xueii once the numerical 
coefficient is changed in correspondence with the different T. For this comparison we 
must quote that CLOUDY ioi z = A gives for sQSO, sQSO+star, sQSO+sDDM and 
sQSO+wDDM models R ~ 17, 37, 160 and 47 kpc respectively with T4 ~ 4.5. 

CLOUDY has produced rather large temperatures for LC which might result 
from thermal evolution of collapsed systems (Miralda-Escude & Rees 1994) and that 
correspond to Doppler broadening in the range 25 km s^^ ^ 6 ^ 32 km s^^, which is 
marginally consistent with the measured values obtained from high resolution QSO 
spectra. This would require gravitational confinement of LC in a colder IGM. 

Considering only the photoionization equilibrium for Hell, it is also possible to 
check the consistency of predictions for different thin absorbers through 
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Xucii II 



(10) 



1 — XHell J Hell CTHell 



Eq. 10 is one of the ingredients which lead to Eq. 9 and involves only the density and 
temperature of the medium so that it is valid also for homogeneous IGM. 

Since for LLS and LC CLOUDY gave the results of Tables 1 and 2, it is clear that 
the discrepancies in the use of Eq. 10 to compare these absorbers must come from 
self-shielding effects in LLS. Thinking in principle that UV flux and ionization cross 
section are common properties, the ratio of Eq. 10 for LLS and LC should be 

[XHeIl/XHeIIl]LLS [XHelll I XHeIl\LC = i^Hi^LS / ^H^c) {Tlc/TlLS?^^ ■ (H) 

If we take the values of Tables 1 and 2 for QSO, QSO+star and QSO+DDM 
models at 2; = 2, 4 being typically T4 ~ 2 for LLS and 5 for LC, the left-hand side of 
Eq. 11 is larger than the right-hand side in an amount which is explained by the same 
reduction of JrcII for LLS that we discussed in connection with the discrepancies of 



Therefore we conclude that the approximate expressions based on photoionization 
equilibrium work well for thin absorbers so that we may use them to compare LC 
with the less dense homogeneous IGM. 

4. Homogeneous Intergalactic Medium 

We anticipate that the results of this Section cannot be too precise due to the 
large uncertainties in the properties of homogeneous IGM. 

As a first step we consider GP effect for HI which, if the ionization is due entirely 
to UV radiation, will correspond to the optical depth 



where VLigm is the fraction of critical density in the homogeneous IGM. 

^From the bound (Giallongo et al. 1994) aX. z = 4.3 r^f < 0.02 we will obtain 
the bounds on density of IGM for different UV models which then, through GP for 
Hell, will allow to estimate XhcII ■ 

We start with the sQSO+star model with Sl = 100 at 2; = 4, taking T4 ~ 2 
which corresponds to inhomogeneous photoionization (Miralda-Escude & Rees 1994). 
We adopt the value of the Hubble constant h = 0.7 which seems to emerge recently 
(Preedman et al. 1994; Riess et al. 1995; Tanvir et al. 1995), though also smaller 
values have been obtained (Tammann et al. 1996; Branch et al. 1996), that suggests 
the presence of the cosmological constant (Krauss & Turner 1995). Since for 

~ 0.6 plus CDM, model denoted as ACDM, the ratio H{A)/Hq ~ 7.1 is obtained 
(Reisenegger & Miralda-Escude 1995), the estimation flicM < follows from the 
bound of Eq. 12 where we have taken the average of the quoted values of the baryonic 
density Qb — 0.015//i^. It is interesting that this agrees with simulations for ACDM 
giving 80% of baryons in collapsed form (Miralda-Escude et al. 1995). 

We now pass to GP for Hell whose optical depth is 



Eq. 9. 




(12) 



Q 




0.015 )\hj\ 4 J H{z) 



VHell 



(13) 



where yneii — nHeii/f^H- The observations (Jakobsen et al. 1994; Tytler et al. 1995; 
Davidscn ct al. 1996) and their interpretation (Shapiro 1995) allow to estimate 
T^fi ~ 1 at 2; = 3.3 . 

With the above parameters and satmating the bound for VLigm the QSO+star 
model gives from Eq. 13 yneii — 1-7 x 10~^ and, using the cosmological ratio between 
H and He, XhcII — 2 x 10~^. A similar prediction results from the addition of modified 
decaying neutrinos to QSOs. For QSOs alone the flux Jri -22 ~ 2 at 2; ~ 4 would 
produce instead a larger fraction XhcII — 3 x 10~^. 

For the standard QSO+sDDM alternative if one takes the averaged flux 
J HI -22 ~ 30 at 2; = 4, T4 ~ 1 due to the lower temperature of homogeneous photoion- 
ization and the rest of parameters as required by the decaying neutrino model (Sciama 
1990a) n = 1, H[z)/ Hq = (1 + 2;)'^/^ , h = 0.56 (which would coincide with the mea- 
surement of Tammann et al. 1996 and Branch et al. 1996), from Eq. 12 it turns out 
^iGM < l^s- The saturation of this bound for fljoM would be in agreement with 
the delayed formation of structures predicted by HDM which requires fi/cM > 
(WiUiger et al. 1994). In this way one would obtain from Eq. 13 yHeii — 0.9 x 10~^ 
and XHeii ^ 1.1 X lO'^. 

We now compare these results for IGM with those for LC of Table 2 using Eq. 10 
from which 



Taking as average the properties of IGM at 2; ~ 4 it seems clear that the spectrum 
of QSOs alone is excluded because, since nHjQM ~ 10~^ Qjgm /^b cra~^ (Madau & 
Meiksin 1994), inserting the rest of parameters in Eq. 14 the left-hand side is one 
order of magnitude smaller than the right-hand side. This conclusion is consistent 
with the fact that the UV spectrum must be softer than that of QSOs to agree with 
the estimations of GP as seen in Eq. 1. 

With the QSO+star model the left-hand side of Eq. 14 increases because XhcIIigm 
is smaller, and the right-hand side decreases because Uhjom larger so that the 
disagreement diminishes. A similar situation applies to the addition of modified 
decaying neutrinos to QSOs. One must note that this conclusion is independent of 
the chosen value of hH(z)/Ho since both XHeii and uh for IGM turn out to depend 
on the square root of it. 

For the sQSO+sDDM alternative the comparison of both sides of Eq. 14 is reverted 
because, to the further modifications of XhcII and hh for IGM, one must add its 
predicted smaller temperature so that the left-hand side becomes almost twice larger 
than the right-hand side. This indicates that Sl may be too large in this case. 

It is interesting to note that for the QSO+star model the increase of Sl should 
improve the agreement of Eq. 14. In fact with our reduction of 50% of the flux due to 
QSOs which gives Sl = 200 at 2; = 4, XHeii for IGM increases in 100% because one 
would expect r^^jj ~ 2. But looking at Eq. 9 XHeii for LC increases more because, 
apart from the influence of JhcII, a smaller global flux decreases R and T giving a 




(14) 



in 



further enhancing factor. However the improvement is not fast for Sl = 200 as seen 
in Table 2. 

We must note that the same procedure of using a weak QSO contribution cannot 
be apphed to the modified decaying neutrino because, in doing so, one would obtain a 
ratio Jhi/ JhcI at 2; = 2 which would clearly exceed the bound of Eq. 3. Therefore the 
only way to improve the comparison of the QSO+wDDM model with homogeneous 
IGM would be to accept a decline of QSO contribution for 2; > 3 more marked that 
the generally estimated exp[—Q.^%{z — 3)] assumed here. 

Apart from the modification of decaying neutrinos obtained diminishing slightly 
its mass, one could keep this value to allow Nl ionization and increase its lifetime. 
Relaxing the condition that neutrinos ionize alone the HI in the Milky Way (Dodelson 
& Jubas 1994) one might assume the v to be more stable in one order of magnitude 



Several models for massive decaying neutrinos give the dependence of lifetime on 
mass and magnetic moment 



and in particular the minimal supersymmetric standard model with couplings which 
violate the R parity allows easily (Roulet & Tommasini 1991) m,^ ~ 29 eV and 
H ~ 10~^^/iB . To increase the lifetime to Ty ~ 10^^ s one must decrease slightly these 
coupling constants and, to keep the value of niy , correspondingly increase the scale 
of supersymmetric partners above 100 GeV. This is perfectly possible in the range of 
the theoretically admissible parameters. 

The averaged DDM flux would be similar to that of the wDDM previously de- 
scribed so that again, added to that of sQSO model, the properties of LLS would be 
reproduced but the consistence between LC and IGM would be difficult because of 
not high enough Si- 
ll seems that the most convenient shape of UV spectrum for a general agreement 
is one where the flux decreases gently between HI and Hel ionization frequencies 
and then drops abruptly towards the Hell edge. A modification of DDM that might 
give this result is to further increase Tj, above 10^^ s and simultaneously double the 
neutrino mass to allow the ionization of both HI and Hel, which would be admissible 
\l h > 0.74 (Bradford & Hogan 1996), but since this imphes a major change of the 
model we will not analyze in detail its consequences here. 

Obviously all our discussion related to homogeneous IGM is based on assuming 
photoionization as seen in Eq. 12. Whereas for densities n ~ 10^^ — 10^^ cm~^ 
coUisional ionization is not important (Haardt & Madau 1995), it may be relevant in 
low density IGM if the temperature is higher than what assumed here. In this case 
the conclusions regarding UV sources would not apply to IGM. 

5. Conclusions 

Using a photoionization code for absorbers we have seen that the QSO+star model 
for UV sources is able to reproduce the bulk of properties of homogeneous IGM, LC 
and LLS, requiring that S"!, > 100 at 2; ~ 4 in agreement with what emerges from 
recently observed mctallicity abundance. For the comparison with IGM it would be 
crucial to determine its temperature with more precision. These same properties are 



I.e. Ti, ~ 
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not easily reconciled assuming the addition to QSOs of neutrinos of mass ~ 29 eV 
which would complete the critical Universe mass with only HDM. In fact if one takes 
the lifetime ~ 2 x 10^^ s, there is a rough agreement between LC and IGM but 
properties of LLS are not reproduced. On the other hand if one reduces the neutrino 
mass to 27.6 eV or increases its lifetime in one order of magnitude, the LLS difficulty 
is solved but the matching of LC and IGM properties becomes questionable unless 
there is a very fast decline of QSOs for z > 3. 

It is interesting that the original model of decaying neutrinos gives, through GP 
bounds, the large density of non-collapsed baryonic matter predicted by HDM. In 
the same way the QSO+star model is consistent with the scenario of CDM with 
cosmological constant Qa = 0.6 convenient to explain the X-ray emitting intracluster 
gas, as well as old galaxies at large redshift (Krauss 1996), but in difficulty with 
recent measurement of the deceleration parameter (Dodelson et al. 1996). Therefore 
it seems that when GP effect and properties of absorbers will be better determined, the 
observed properties of the reionization of Universe will give definite hints, regarding 
the models which originated the structures, at a redshift intermediate between the 
one of recombination age and that of formation of bulk of galaxies. 

We are deeply indebted to Prof. D.W. Sciama for comments related to modified 
decaying neutrinos. S.S. thanks heartily the hospitality at the Osservatorio Astro- 
nomico di Roma where part of this research was performed. L.M. acknowledges partial 
financial support from CONICET of Argentina through grant PID 339650092. 
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